Feedback stabilization of pure states in quantum transport 
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We propose a feedback control scheme for generating and stabilizing pure states of transport 
devices, such as charge qubits, under non-equilibrium conditions. The purification of the device 
state is conditioned on single electron jumps and leaves a clear signal in the full counting statistics 
which can be used to optimize control parameters. As an example of our control scheme, we are 
presenting the stabilization pure transport states in a double quantum dot setup with the inclusion 
of phonon dephasing. 

PACS numbers: 73.23.Hk, 73.63.Kv, 85.35.Gv, 03.65.Yz 



I. INTRODUCTION 



The use of feedback to stabilize desirable structures or 
behaviors that would otherwise be unobservable is a topic 
of considerable interest in classical dynamics^. Feedback 
concepts have been extended into the quantum realnt^ii^, 
and have been shown to be effective in e.g. prolonging 
the life of Schrodinger cat states^, increasing the rate 
of state reduction under measurement^, and for state 
purification^. Feedback schemes have also been proposed 
in the solid state, both theoretically^ and in a recent 
experimenl)^. In these latter examples, a quantum trans- 
port device (e.g. quantum point contact, single-electron 
transistor) acts as a detector of the system (typically a 
qubit) which is to be controlled. By feeding information 
gained from the measurement device into the system in 
a control loop, one can purify the quantum state of the 
qubit. 

In this work, we take a different approach and couple 
information gathered about the charge transport through 
a device back into the device itself. More specifically, 
we monitor the flow of electrons from the source lead 
to a device in the Coulomb blockade regime, as in full 
counting statistics (FCS)^"— , and conditioned on elec- 
tron jumps, perform unitary operations on the device. 
Using this scheme, we show how such feedback can lead 
to a purification of the device wave function, despite it 
being in non-equilibrium (this contrasts with the situ- 
ation without control, where the state of the system is 
decidedly mixed, as one might expect). We view this be- 
havior as a genuine state stabilization as these pure states 
are implicit in the behavior of the device without feed- 
back; but only with the feedback do they become stable. 
Furthermore, we show how successful purification of the 
transport state also leads to a distinctive signature in the 
statistics of electron flow through the device. We describe 
how this information may be used in the feedback loop 
to optimize the stabilization. After presenting a general 
theory of state stabilization of non-equilibrium quantum 
transport, we consider the example of a double quantum 
dot (DQD) (Fig. [T]) and show that the stabilization of 
coherent delocalized states is possible. 




FIG. 1. Main part: After an electron jumps into the double 
quantum dot (DQD) (measured e.g. with a quantum point 
contact QPC) a pulsed control operation modifies e.g. the 
gate voltages Va and rotates the electron into a different 
state. In an effective model this control operation modifies 
the left jump operator. The DQD is described by e—eL—en 
the detuning and Tc the coupling strength between the dots. 
The electrons tunnel into/out the left /right dot with the rate 
r^^/fl. Inset: For sufficient control operations the DQD is 
stabilized. The transport can then be described by a single 
resonate level model with an effective right tunnel rate 



II. GENERAL THEORY OF THE CONTROL 
SCHEME 



A. Transport model 

Wc start by considering a general transport Hamilto- 
nian, H = Hs + -I- iJx, composed of system (e.g. 
quantum dots (QD)s), leads and tunnel-coupling parts. 
We consider two leads with non-interacting Hamiltonian 



H 



Cak 



with 



the annihilation oper- 
ator of an electron in lead a ~ L,R (left, right) with 
quantum numbers k. We assume strong Coulomb block- 
ade and a bias configuration such that at most one ex- 
cess electron can be in the system at any given time. 
The state-space of the system is thus spanned by the 



2 



empty state |0) and N single electron states typ- 
ically not eigenstates of system Hamiltonian Hg. Fi- 
nally we assume that the leads are each coupled to one 
and only one localized system level (|L) and \R) respec- 
tively) such that the tunnel Hamiltonian reads Ht = 
EktRkclkDB. + tLkci,D[ + H.c. with Dl=\L){0\ and 

DlM0){Rl 

In the high-bias limit with tunnel rates Fq, — 
^'''Sfc I^QA;P'5(e — Cafc) assumed constant, the behavior 
of the system density matrix can be described with a 
Markovian master equation of Lindblad forroi^: 



(1) 



P = Wp= (Wo + Jl + Jr]p 



with jump super-operators JaP^^aD^pDa, which trans- 
fer an electron to/from lead a = R/L, and free Liou- 
villian Wo, which describes the evolution of the system 
without electron transfer. This latter assumes the form 
Wqp ~ —i — pH^^ (we set h = 1), where 



(2) 



is an effective non-Hermitian 'Hamiltonian' operator for 
the system. 



B. Electron counting and control 

The evolution of the density matrix under Eq. ([IJ can 
be interpreted in terms of trajectories^^. The density 
matrix at time t is then given as a sum of terms such as 



t rt2 

dt2 / dma{t~t2)jRna{t2-ti)JMti)p{t^), 
Jo 



(3) 



which describe all trajectories in which an electron jumps 
into the system at time ti (iTl), out again at time t2 
with evolution in between described by the no- 
jump propagator r2o(t) = e^"*. As in FCS experiments, 
we assume that our measurement device is sensitive to 
single electron jumps. Conditioned then on the detection 
of an electron jump from the left lead, we set the control 
electronics to implement a control operation on the sys- 
tem which we consider here to act as an instantaneous 
unitary operation on the system^. In an experimental 
setup for the transport through QDs these control oper- 
ations could be achieved, for example, with pulsed gate 
voltages. This scheme modifies the quantum trajectories 
such that each left jump operator is followed by a con- 
trol operation; we thus replace Jl with = ClJl in 
all trajectories like Eq. ([3]). Re-summing all the trajecto- 
ries, we arrive at an effective Liouvillian for the evolution 
of the density matrix with control that is the same as W 
in Eq. ^ but with Jl replaced by jf^ . 



C. Feedback Stabilization 

We now describe how this control schema can be used 
to stabilize certain states. The effective "Hamiltonian", 
Eq. ([2), has right and left eigenstates H\tpj) = £j|V'i) 
and {i^jlH ~ Sj{'il'j\, which, in general, are non-adjoint: 

{tpj\ 7^ i\fpj))^ , and can be used to write H in its diagonal 
basis as 



N 

E 



(4) 



where we have separated the empty state since it is not 
coupled to other states by H. The eigen-operators of the 
free Liouvillian are Pjk=\ipj){4'k\ such that WoPjk= — 
i (Ej — e^) pjk- The diagonal matrices pjj represent pure 
states and obey yVoPjj=2'^{£j)pjj ■ The empty state is 
one such eigen-operator with Wo Poo = ^TlPoo- The 
remaining pure density matrices pjj; j > are states 
that we will seek to stabilize with our feedback scheme. 

To effect the state stabilization, wc choose the control 
operation Cl to rotate state \L) into the desired eigen- 
state of H. The effective left jump operator therefore acts 
as Jl Poo ~ ^LPjj, where the F^ forefactor arises from 
the normalization and the fact that the empty state de- 
cays with rate F^ with or without control. Once in state 
Pjj, no transitions out of this state are induced by Wo- 
Rather, the only thing that happens to an electron in this 
state is that it leaves at a rate 'y^'^ ~ — 23(ej). By jump- 
ing directly into one of the eigenstates of the free Liouvil- 
lian Itpj) the dynamic of the system is then determined 
solely by the vacuum \ipo) and the other states are 
decoupled. In the basis {poo^Pjj} the Liouvillian of the 
feedback controlled system in the above parametrization 
can be written as 



U) 
c 




(5) 



This Liouvillian corresponds to an effective single res- 
onant level, with the effective tunnel rate 7^' and the 
stationary state 

1 



Pstat 



0) 



Fl 



7fl 



(6) 



In the limit F^ ^ 7)5; the stationary state reduces to 



lim 



(7) 



and the system is thus stabilized in the pure state \ijjj). 

In general then, a system with N internal states has N 
pure states that can be prepared in this manner. These 
states depend on the internal system parameters as well 
as on Fjj. In the limit of F^— >0, the stabilized states 
will be eigenstates of Hs- For finite Tji, however, these 
states are, in some sense, non-equilibrium pure states, 
since they depend on the rate F^j and can be quite dif- 
ferent from the eigenstates of the bare Hamiltonian Hg. 
Needless to say, the states are unstable without control. 
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III. FEEDBACK STABILIZATION OF A 
DOUBLE QUANTUM DOT 



In the following, 



we discuss the N= 
1 , 



■2 example of a 



charge qubit with Hamiltonian Hs = ^^ctz + Tc(Jx, with 
a, = |L) {R\, {R\ + \R) {L\, t^eL-^R the dc- 

tuning and Tq the coupling strength between the dots. 
In the infinite bias limi t^'^i^^i"'^^ the transport through 
the DQD is described by Eq. ^ with d[=|L)(0| and 



D]j=|0)(i?|. The corresponding Liouvillian of the DQD 
is shown in appendix 1X1 



A. Zero detuning e = 



First we consider the case without detuning e=0 for 
simplicity. The right-eigenvalues of iJ are eo=— with 
eigenvector |-0o) = |O) and 



16T2 



(8) 



with the eigenvectors j?/'^) = M^^\(iVii =F k)|L) + 

4rc|i?)] , normalized with M^^^/wF^^TWr^W ■ 

In the following, we consider the simplification ^ 
FfljTp which allows us to effectively project out the 
empty dot state |i/'o) and to assess the purity of the 
transport qubit in the usual Bloch representation with 
{ax}=2^{pLR), {ay)=2'^{pLR) and {iJz)=Pll- Prr- 

The pure-state density matrices pipip = |'(/'=F)(V'Tl to be 
stabilized correspond to two of the occupied eigenstates 
of the free Liouvillian Wq- Fig. [^Ja) depicts these stabil- 
isable pure states on the surface of the Bloch sphere. For 
Fi^>4Tc they describe partly localized states with 



i R. 



and 



AT, 



c 



(9) 



and for Tii<ATc completely delocalized states with 



4T( 



and 



c 



4rc' 



(10) 



with K=^F^ - ler^. At Tr=ATc the eigenvalues Ezp 
are degenerate and the corresponding eigenstates cross 
each other with {az)=G, {ax)=0 and ((Tj,)=1. 

Without control an electron simply tunnels into the 
left dot (and thus into a localized, pure state). The 
goal of the control operation is to rotate this state 
into one of the pure states p±±=\ip±}{'>p±\- In gen- 
eral, we can consider every qubit rotation as a possi- 
ble control operation Ci,=e*^"'°'^ where S corresponds 
to the Pauli matrices in Liouville space 'Ep~—i2[(T, p] 
with the usual Pauli matrix vector cr, the unit vec- 
tor 'nQ~{nri:,ny,nz}={sin9,0,cos9} that determines the 
angle of the rotation in the (T^-crz plane, and the 
control strength parameter 9c- A straightforward 
calculation (see appendix [B]) leads to the control 
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FIG. 2. (a) Stabilisable double quantum dot states pqiip 
on the surface of the Bloch sphere, and stationary state pwo 
without control, for tunnel coupling Tc £ {0, oo} . The tri- 
angle marks the state which is stabilized in (b). (b) Time 
evolution of Bloch vector length | (cr) | with control operation 
switched on at t = SF^^ for different temperatures T. Param- 
eters: Tc ~ 2r_R, phonon coupling g — 8- 10"'', Debye-cutoff 
ujc = 500r_R and control parameters 6 = arccos(3i/7/127) 
and 6c = arcsec(8-\/2). 



parameters 9 
for Ffl>4rc, 



-tt/2 and 
and 



ic~ arccos 



arccos 



(V(rfl±K)/2Ffl 



-^2 
R 



32r2 



and 



Wc^arcscc (4V2Tc/Fi?) for Tr<ATc, corresponding to 
the stabihzation of p±±. The upper left plot of Fig. [3] 
shows these control parameter on the control plane. 

The fidelity of our control scheme is characterized by 
the length of the Bloch vector | (cr) | , obtained by solving 
Eq. ([U for the DQD model. Fig. [IJb) (sohd hue) demon- 
strates the effect of the control operation Cl on the DQD 
Bloch vector length | (cr) | — switching on the control (at 
time t = 8F^^ in the plot) leads to a rapid purification 
of the transport qubit state, even though electrons con- 
tinuously tunnel through the DQD. 



B. Finite detuning e 7^ 

With finite detuning it is possible to stabilize all states 
on the Bloch sphere which have a positive {<7y) compo- 
nent. The exact expressions for the stabilizable eigen- 
states with e 7^ can be found in appendix [Bl 

Fig. [3] (right) shows the projection of the stabilizable 
eigenstates on the {(yz)-{cx)- plane of the Bloch sphere for 
different detunings for Tc G {0, 00}. Without detuning 
(e = 0) the stabilizable states lie on the axis of the plane 
which can be also seen in Fig. [IJa). Since all stabilized 
states have a Bloch vector with length one and a positive 
(tTj,)- component, the size of (ctj,) is maximal when the 
states in the {a'z)-{o'x)- plane have the shortest distance 
to the origin. For increasing detuning these states move 
asymptotically towards the unit circle, meaning the max- 
imal size of the (ay) component decreases for increasing 
|e|. Each of the four quadrants of the plane is covered by 
a different solution: The I- quadrant by the solution 

for e > 0, the II- quadrant by the p solution for e < 0, 

the III- quadrant by the solution for e > and the 



4 




FIG. 3. (left) The control operations needed to stabilize states 
at different detunings e for Tc G {0,00} are shown on the 
control plane. The control strength 9c is determined by the 
length of a vector between a point on the control plane and 
the origin. The angle between the positive part of the Ocriz- 
axis and the vector is the control angle 6. (right) The cor- 
responding stabilizable states are shown as projection on the 
{a,)-{ax}- plane of the Bloch sphere. 



IV- quadrant by the p++ solution for e < 0. The con- 
trol operation necessary to stabilize the states for each 
detuning are shown on the left side of Fig. [3l 



C. Phonons 

In any realistic experimental setup, our ideal coherent 
feedback control scheme will be disturbed by the cou- 
pling to a dissipative environment. In order to make 
quantitative predictions, we therefore introduce an ad- 
ditional coupling of the DQD to a bath of thermal 
phononsil with wave vector Q and energy ujq. We de- 
scribe the electron-phonon coupling via a (pscudo) spin- 
boson Hamiltonian^*' 



phon 



UJQaQQQ 



(11) 



that leads to a modification of the free Liouvillian Wo 
in our master equation Eq. ([Ij. These modifications are 
included in the free Liouvillian shown in appendix [A] We 
assume that an electron in the DQD couples to a plionon 
bath with an ohmic spectral density with exponential 
cut-off at Debye frequency Wc- Without detuning (e = 0) 
the corresponding dephasing rates at temperature T (we 
set the Boltzmann constant fc^ = 1) are 



7p = 7 coth 



Tc 
T ' 



7 = gnTce' 



76 = 0, (12) 



where 5 ^ 1 is a dimensionless coupling parameter. 

Fig. IHb) also shows the evolution of the fidelity with 
the inclusion of the phonon bath (dashed and dotted 
line) at finite temperature. The lengths of the stabi- 
lized steady-state Bloch-vectors with phonons is below 
one I (cr) | < 1 since the phonon dephasing limits the sta- 
bilization of the pure states. 

Fig. SJa) shows the fidelity of the steady state for e = 
as a function of tunnel coupling Tc- The stationary state 
without control {9c=0) is always mixed (except for triv- 
ially at Tc = where |(<t)|=1 as then an electron tunnels 
into the left dot and never leaves it), with | (cr) |— >-0 at very 
large Tc- In contrast the control-stabilized states p±± 
remain pure in absence of phonons (|(<t)|=1 for g = 0). 
These pure states are differently affected by the phonon 
dephasing. This striking behaviour of the two phonon 
branches for the control-stabilized states in Fig.|4l^a) 
can be traced back to the character of the eigenvalues 
e^, Eq. ([8]), and the corresponding cigenstatcs of 
H. Starting from large Tc ^ the Eip have a much 

larger real than imaginary part and {tp-) (|V'-i-)) essen- 
tially corresponds to the ground (excited) state of the 
isolated DQD. In this regime of large level splitting, spon- 
taneous emission of phonons from \ip+) to \4'~) dominates 
phonon absorption, with the purity of p++ strongly re- 
duced as compared with the one of p--, which is shown 
by the explicit expressions for the fidelities |((T-|-)| of the 
two states, 




1 ± tanh- 



Tc 
T 



(13) 



In particular, the shape of the p++ branch for large Tc 
and at low temperatures is entirely determined by the 
spectral density of the phonons oc Tce~^"^'^/'^= and can 
therefore be used as a spectroscopic tool to determine the 
dephasing rate 7p, Eq. ((T^ . 

Moving to smaller tunnel coupling Tc, the two broad- 
ened levels of the DQD start to overlap at Tc=Tii/4: at 
which point the e^, Eq. ([8]), become purely imaginary, 
and the two merged branches split again at still smaller 
Tc- In absence of phonons, our control scheme coherently 
rotates the DQD into the pure states \ip±), i.e., |((t±)|=1 
for g=0. The presence of phonons ((?>0), however, leads 
to an un-controUcd decay of 1-0^) that can be described 
by — )■ Tii+2"fp in the imaginary eigenvalues Eq. ([5]), 
which for Tc^Tr then become e+wO (very slow decay) 
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FIG. 4. Loss of stabilization of the control stabilized pure 
states pzpip due to phonons illustrated by the length of the 
Bloch vector as a function of Tc for different detunings (a) 
e = (without detuning), (b) e = ±^r_R and (c) e = ±10r_R. 
Without control {9c = 0) the length of the Bloch vector 
decays to zero for large tunnel couplings Tc- Parameters: 
Vl — > oo, phonon coupling g — 4,- 10~*, luc/Tr — 500, 
fcsT = 50rfl. 
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FIG. 5. Fano factors F^"' for the double quantum dot with 
zero phonon coupling as a function of the control strength 
9c at fixed qubit rotation angle 9 = ^. Stabilization of a 
pure state occurs at the value of 9c where F'"-* = l;Vn > 0, 
corresponding to Poissonian statistics in the limit of 3> 
TcTr considered here, (e = 0) 



large Tc — > oo close to {<Tx) = ±1. The phonon dephas- 
ing has the strongest effect on the states which either 
lie close to (a^) = —1 or on the states with (ax) ^ 1- 
The control operation in the first case rotates the elec- 
tron that originally tunnels in the left dot into a state 
which is almost entirely localized in the right dot and in 
the second case the electron is rotated into a state which 
corresponds to the excited state of the isolated system 
(see discussion for e = 0). In both cases the states are 
strongly effected by dephasing. In the first case this ef- 
fect vanishes for decreasing temperatures but in the sec- 
ond case spontaneous emission even at zero temperature 
destabilizes the states, with the exception of states above 
the Debye cut-off (these states are only weakly effected 
by the phonon dephasing). This knowledge explains the 
anti-crossing and crossing in Fig.|4ljb) and (c). For e < 
the solution stronger destabilized by phonons is the same 
{p++) for small and for large Tc which leads to an anti 
crossing and for e > the solution stronger afi'ected by 

dephasing switches form p for small Tc to /?++ for 

large Tc resulting in a crossing. 



and — i{rji/2 + jp) (fast decay). This 'repulsion' of 
imaginary levels is analogous to the Dicke spectral line 
effecti^ii^, i.e. a splitting into a sub- and a superradiant 
decay channel, and is again confirmed by the fidelities 

|(<T+)|«1, (14) 

valid for Tc't:TR. 

The fidelity with finite detuning as a function of Tc is 
shown in Fig.lH^b) and (c). The main features of the plots 
are the same as in the case without detuning Fig. lU^a). 
However if the detuning between the dots is positive e > 

the P++- and the p - branch cross each other, while for 

negative detunings e < they preform an anti-crossing. 
In general the stabilizable states (see Fig. El^right)) for 
small couplings Tc — !• lie close to (az) = ±1 and for 



D. Full counting statistics 

Collecting all trajectories with n left-jumps as partial 
density matrix p^"-* , Eq. ^ implies the number- resolved 
master equation for the system with control 

= (Wo + Jb.)p^"^ + JEp^"''^- (15) 

With the Fourier transform p{x) = e'"-^p^"-' we ob- 
tain the ^-resolved master equation 

Kx) = (Wo + Jr + J-f e*'^)p(x)- (16) 
Using the current cumulant generating functioniSi"— 

^(x,t) = ln(Trz5{e(^''+-^«+'^"^'^)(*-*°)p(to)}), (17) 
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where Tr£){- • • } corresponds to the trace of density ma- 
trix, we calculate the 71-th order zero frequency current 
correlation 



d r)" 

We define the zero-frequency Fano factors as 



(18) 



(19) 



the n-th order current correlation functions, normalized 
by the stationary currenlii^ {{^'•^'') = {!))■ 

Fig. [5]shows the zero-frequency Fano factors F^") with 
n = 2, 3, 4 for our DQD example as a function of control 
strength 9c with fixed static dot parameters and control 
angle ^ = f (again in the limit F^ — >■ 00). Whilst the 
cumulants exhibit a rather complex dependence on 6c 
at exactly the point where the pure qubit state is stabi- 
lized, all Fano factors are simply = 1. This can be 
understood by recalling that at the stabilization point, 
the system is equivalent to a single resonant level in its 
transport properties, cf. Eq. ([5]). The effective outgoing 
rate for our DQD example for e = reads 



i(rfl±yrf" 



rfl,<4rc 

16r2) Tr>4:Tc 



(20) 



The cumulant generating function of this effective model 
in the long time limit is wcll-knowoi^^ 



-7i'V V^(rr^^jf)^74rl^^) . 

(21) 



Which in the limit F^ — > 00 reduces to the cumulant 
generating function of a Poissonian process for which all 
Fano factors are equal to unity. 

This information that the stabilized state has such a 
distinctive FCS can be used to locate the control op- 
eration required for stabilization even in the absence of 
knowledge of the precise model parameters. Knowledge 
of the electron jumps permits the FCS of the transport 
to be calculated. For a given set of, presumably non- 
stabilizing, control parameters, this FCS can be com- 
pared with an ideal Poissonian signal as a target. The 
control parameters can then be adjusted and the FCS 
compared again. This can be implemented in an in situ 
(classical) closed feedback loop to optimize the control 
operation such as to permit stabilization. With elec- 
tron jumps occurring, e.g., on time scales of milli-seconds 
in state-of-the-art FCS experiments^^, no ultra-fast elec- 
tronics is necessary for information processing. The opti- 
mization should also work if the condition of large F l can 
not be assumed, in which case the target signal would no 
longer be Poissonian but given by the cumulants of the 
two coupled Poisson processes (left and right tunnel bar- 
rier) of the resonant level model Eq. (HI]), Refj^iii^. 



IV. CONCLUSIONS 

We have presented a general theory for the stabiliza- 
tion of pure states in a non-equilibrium electron transport 
setup. The stabilized states have the features of a single 
resonant level including the FCS, which can be used to 
detect these pure states. We demonstrated this stabi- 
lization on the example of a DQD and have shown that 
half of the states on the Bloch sphere can be stabilized. 
Moreover the interaction with an dissipative environment 
due to the coupling to a phonon bath was discussed. Fur- 
ther generalization of the control stabilization to systems 
with many body states and more leads enable to use this 
stabilization scheme for arbitrary Coulomb blockade sys- 
tems. 

We are grateful to G. Kiefilich, G. Schaller, 
A. Beyreuther and K. Mosshammer for useful discus- 
sions. Financial support by DFG projects BR 1528/7, 
BR 1528/8, GRK 1558 and SFB 910 is acknowledged. 



Appendix A: Double quantum dot Liouvillian 

We present the DQD Liouvillian with phonon 
coupling in the pseudo-spin basis with p = 

{POO,"-occ, (CTx), (o-y), (CTz)} = {PO,PLL + PRR,PLR + 

Prl, }{plr - Prl),Pll - Prr}- 
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The left jump operator without control J7l can be calcu- 
lated by setting 6'c' = in J'^ . For an Ohmic spectral 
density with exponential cut-off at Debye frequency tOc, 
the phonon dephasing rates are 
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where g ^ 1 is a dimensionless coupling parameter, 
/3 = (ksT)^^ is the inverse of temperature T times Boltz- 
mann constant fcs and A y/e^ + iT^. In order to find 
the stabihzable states the phonon coupling has to be set 
to zero (g = or 7 = 7p = 7b = 0). 



Appendix B: Stabilizable states 

The DQD- eigenstates we seek to stabilize are p±±. 
The control operation which stabilizes these states are 

jC 

determined by Yf^poo = P±±- In general, the right eigen- 
values of H of the DQD which correspond to occupied 
states are 

e^^^{-iTRT v/4e2 + AieVn - Tj, + 16T^) . (Bl) 



It is possible to calculate the stabilizable states di- 
rectly from the left- and right- eigenstates of H for e 7^ 0, 
this is cumbersome due the imaginary part under the 
square root. In practice it is easier to diagonalize Wg 
directly in the spin basis. Due to normalization the riocc 
components of the diagonalized occupied Wq- eigenstates 
are equal to one. Since we are interested in the two oc- 
cupied eigenstates corresponding to pure states, the nor- 
malized (cTi), i £ {x, y, z} components have to be real and 



\/ {<^x)^ + i'^y)^ + i'^'ij)^ = 1- The eigenvalues of these 
state are A^^ = -^Tb T ^ - 4e2 - 16T'^ + Tw 
with Tw = V(4e^ -f + 16T^)2 - 64r^r^ with 



+ r|, + i6r2 - r 



y/r^ - 462 - lerg + r, 



(B2) 



In the limit e — >■ these eigenstates reduce to Eq. ^ for 
Vr > ATc and to Eq. ^ for < 4Tc. In order to 
calculate the control operations, these components have 

jC 

to be equal to matching components of -^Pqq: 

K) =sin(20)sin(0c)', 
(ay) =sin(0)sin(20c), 
{a,) =cos(6')2 +cos(26lc)sin(6')2. 



(B3) 



Careful analysis shows: For e > one has 25(e^) = A^^, 
meaning the eigenstate pipip corresponds to the eigen- 
value Azpzp of the free Liouvillian. While for e < we 
have 29(e;p) = A±± and the eigenstate corresponds 
to the eigenvalue K±± of the free Liouvillian. 
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